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Abstract 

The thermally driven confinement- deconfinement transition exhibited by lattice quantum 
electrodynamics in two space dimensions is re-examined in the context of the statistical 
gauge-fields common to anyon superconductors and to spin-liquids. Particle- hole excita- 
tions in both systems are bound by a confining string at temperatures below the transition 
temperature T^,. We argue that coincides with the actual critical temperature for anyon 
superconductivity. The corresponding specific-heat contribution, however, shows a smooth 
peak just below Tc characteristic of certain high-temperature superconductors. 
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The discovery of high-temperature superconductivity-^ has generated a number of 
new theoretical approaches to the problem of strongly interacting Fermi systems in two 
dimensions. Among these, the proposals of spin-liquids in two-dimensional (2D) Mott 
insulators and of anyon superconductivity in doped 2D Mott insulators are perhaps the 
most novel. The former state is characterized by short-range spin correlations while the 
latter is essentially a 2D superconductor that generally breaks time-reversal and parity 
symmetries."^ In addition, both the spin-liquid state and the anyon superconductor have a 
gap for Fermi excitations.^ 

Parallel technical developments have evolved in the arena of ^7(1) gauge-field de- 
scriptions of strongly correlated electron systems in two dimensions. In particular, the 
2D spin-1/2 antiferromagnet has been treated within both the Schwinger-boson scheme^ 
and the corresponding long-wavelength description given by the CPi model, ^'^ where the 
gauge-field that appears in such case measures chiral spin-fluctuations.^ In addition, the 
random-phase approximation (RPA) results for anyon superconductors^ have been recov- 
ered by Chern-Simons (CS) gauge-theories in the mean-field approximation,^'^ where a 
statistical gauge-field is introduced to describe the fiux-tube attached to each fermion by 
the CS term. Very similar gauge-field descriptions also exist for anyon superconductors in 
the context of doped 2D antiferromagnets [t — J model). ^^'^^ In both the case of the spin- 
liquid (quantum disordered) phase of the 2D spin-1/2 antiferromagnet and of the anyon 
superconductor, fluctuations in the (statistical) gauge-field are described by an action for 
vacuum quantum electrodynamics in 2-|-l dimensions (QEDs).^'^'-*^^ This does not occur 
in the case of either the 2D Neel phase^ or of the "strange" metal phase of the 2T> t — J 
model, where spin excitations are gapless. 

In this paper, we study the thermal properties of such statistical gauge-fields com- 
mon to anyon superconductors and to spin-liquids on the square lattice. Specifically, we 
re-examine the confinement-deconfinement (CD) transition at non-zero temperature expe- 
rienced by particlc-hole excitations^^'^^ about such groundstates using an effective lattice 
(compact) QED3 lagrangian for the statistical gauge-field. The CD transition is found to 
be dual to the 2D Coulomb-gas (CG) transition in the weak-coupling regime, like in the 
strong-coupling regime. ^^'^^ In fact, both the string-tension and the inverse of the confine- 
ment length-scale vanish exponentially as the transition temperature, Tc, is approached 
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from low-temperature at both the strong and weak-couphng hmits. The specific-heat is 
also computed, where it is characterized by a smooth bump anomaly below reminiscent 
of certain high-temperature superconductors. ""^^ Last, we argue that Tc coincides with the 
actual transition temperature for anyon superconductivity. 

The action for compact QED3 at temperature T is defined by 



where av{x) represents the statistical gauge- field, denotes the lattice difference opera- 
tor, and where x = {xi,X2,xs) is the three-vector that spans the cubic-lattice space-time 
with lattice constant a. Quantum statistical mechanics requires that the gauge-field be 
periodic in the time-like direction^^ xs = icot, where Cq is the zero-sound speed for anyon 
superconductors and the spin-wave velocity for spin-liquids; i.e., a^{x3 + Pa) = a^^^x^) 
mod 27r, where j3 = hu)o/kBT is presumed to be a positive integer, with uq = coa~^ as 
the Debye frequency-scale of the theory. Polyakov has shown that compact QED3 is con- 
fining at zero-temperature in the weak-coupling regime, go <^ 1. Below, we extend 
his instanton-gas analysis to the present case of non-zero temperature, and find that a 
CD-transition occurs at ksTc ~ gohuQ. 

Weak- coupling. In the limit (70 ^ 1, the Villain substitution for the partition function, 
Z = J T>a^e~^ ^ is then valid, yielding the factorization Z = .^gauss-^inst, where Zgauss 
represents the gaussian ( "spin- wave" ) approximation to the action (1), and where the 
corresponding instanton contribution is 



Above, n{x) is an integer- field on the dual cubic-lattice, such that n^i,{x) = en,^xAxn{x) 
is the antisymmetric integer-field dual to the minima of the periodic energy functional 
(1) at a given plaquette {p, u). [The link associated with Axn{x) passes perpendicularly 
through the plaquette associated with n^i^{x).] Again, the field n{x) must be periodic in 
the time-like direction with period Pa; i.e., the sum over time X3 in Eq. (2) is restrict to a 
slab of thickness Pa. After employing the Poisson summation formula on the configuration 
sum above (2), and then integrating over the continuous field corresponding to n{x), we 




(1) 





3 



arrive at the following instanton-gas ensemble:^° 



Z\nat = XI 
{m{x)} 



1 ' 

^(27r)2 y m{x)G^^\x - x')m{x') 



x,x 



(3) 



where G'('^)(x) = {2n)-'^ Jj^^d'^ke'''-^[2d - 2^J^^^cos{k^a)]-^ is the Greens function for 
the d-dimensional hypercubic lattice, and where m(,x) is the integer charge-field of the 
instanton-gas, which is also periodic in the time-like direction. Consider now two instantons 
of charge m{x) and m{x') centered at space-time points x = (r, ^3) and x' = {f',xs'), 
respectively. Furthermore, suppose that they are widely spaced, such that |r — r '| ^ Pa. 
Then since both configurations have period /?a along the time-like direction, their mutual 
action (3) per elementary time-slice is proportional to f3~^m{r)G^'^\f—f')f3~^m{f'), where 
Q('^)(^f) is the Greens function for the square-lattice and j3~^m{r) is simply the linear charge 
density along time. After multiplying the latter action-density by fi to get the action per 
period, we arrive at the following 2D CG partion function for the periodic instanton gas 
ensemble (3): 



{m{r)} 



*-* f,r' 



(4) 



where the pre-factor of (3^'^ above results from counting all possible instanton configura- 
tions fixed at Nm sites in space. Hence, we obtain a transition at ksTc ~ ksTo = QAAg^TiUQ 
that is dual to the 2D CG transition, ^"^ where Tq denotes the temperature at which the 
chemical potential vanishes in the CG ensemble (4). For temperatures below Tg free in- 
stantons exist, while instanton-anti-instanton pairs are bound above Tg. Note that the 
exponentially small number of free instantons that exist near the dual CG transition^^ 
validates the previous assumption of diluteness. 

To demonstrate that the preceding is in fact a CD-transition, we now compute the 
confinement length-scale (string-tension), and show that it diverges (vanishes) exponen- 
tially as Tc is approached from low-temperature. Adapting Polyakov's zero-temperature 
calculation^^'^^ for the auto-correlation function of the electromagnetic fields in compact 
QED3 to the present case at T reduces to the sine-Gordon reformulation of the dual 
2D CG ensemble (4).^^ A straight-forward generalization of his methods then yields that 
the corresponding static auto-correlation function for the statistical electrodynamic field is 



given by Qq"^ j3{hiik)hj{—k)) = 5ij —kikj/[k'^+^~'^{T)], where the magnetic fields in the spa- 
tial directions j = 1,2 are related to the electric field ei{x) = Fq^i{x) by bj{x) = iejiei{x), 
and where CxC^) denotes the Debye screening- length of the dual CG ensemble (4). Hence, 
the confinement length scale Cx(-^) infinite above Tc in the bound- instanton phase, while 
it diverges below Tc like 

ex(T) = aSexp[^/(l - T/Te)^/2] (5) 

in the free instanton phase, with A and B being non-universal numerical constants. A cor- 
responding adaptation of Polyakov's calculation of the zero-temperature string-tension^°'^^ 
a also yields 

a{T)^glhuJo^-'{T). (6) 

In conclusion, correlations between the statistical gauge field have a finite range below 
Tc indicating confinement, while they are of infinite-range above implying gaussian 
[iT'fj.uix) = 0] electrodynamics consistent with the absence of instantons. This observation 
suggests that the specific-heat contribution of the gauge-fields in the weak-coupling limit of 
compact QED3 follows a gaussian (photon) law above Tc, while the temperature-scale 
below which activated behavior sets in the confined phase is on the order of a{0)a. Here, the 
zero-temperature string tension a{0) is given by (6) with a zero-temperature confinement 
length-scale20'2i of ^^(0) = a(^o/27r)e^°"^*/9o . The former energy-scale is exponentially 
small in the present weak-coupling limit, however, indicating that the specific heat is given 
essentially by the same gaussian T^ law below Tc- 

Strong- coupling. Following Polyakov and Susskind,^^'^^'^^ the Hamiltonian formula- 
tion of compact QED3 (1) reduces to the kinetic energy Hq = |^o"^^'^o (x), in the 
regime ^ro ^ 1) where the statistical electric field operator ei{x) = —igQd / dai{x) satisfies 
Maxwell's Eq., AjCj = 0, in vacuum. (The temporal gauge as = has implicitly been 
chosen.) Because of the compact nature of the vector potential, the field-strength operator 
gQ'^ei{x) has integer eigenvalues ni{x) satisfying Aini{x) = 0. Since the latter constraint 
can be easily solved by taking ni{x) ~ eijAjn{r), where n{r) is an integer-field on the dual 
square-lattice, we arrive at the discrete gaussian (DG) model for surface roughening^^"^^ 
given by 

Z= ew\-l(39'oY.^AMf^A- (7) 

{n(r)} \ f ) 
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This model is dual to the 2D CG ensemble, with a transition temperature ksTc = 
0.73gQ?i(jJo. Notice also that the present DG model (7) is equal to the weak-coupling 
instanton action (2) for the special case of static configurations Asn{x) = 0. Below we 
demonstrate that Tc marks the boundary for the CD transition that also exists in the 
strong-coupling limit. ^^'""^^ 

Let us first determine the temperature dependence of the confinement length-scale for 
compact QED3 in the present strong-coupling limit. Given the identity Qq'^ {ek{x)ei(x')) = 
go^kieij{^i'n{r)Ajn{r')), and the fact that the latter is "proportional" to the height- 
height correlation function of the surface-roughening model (7),^^ we find that the spatial 
Fourier transform of the auto-correlation function for the statistical electric-field is given 

^ki^ijkikj /[k'^ + ^^'^{T)], where denotes the Debye screening 
length of the dual CG ensemble. Hence, we recover precisely the same pole obtained in 
the prior weak-coupling analysis, with a confinement length-scale Cx(^) ^^^^ diverges ex- 
ponentially (5) at Tc, but that vanishes linearly with temperature at zero temperature.^^ 
Note that gQ^/3{e{0) ■ e(0)) jumps from zero to unity upon heating through the transition. 

Confinement properties are also directly probed by the string-tension. To compute 
this quantity, it is convenient to return to the original strong-coupling Hamiltonian that is 
known to be dual to the 2D XY-model. Specifically, Polyakov and Susskind have shown 
that the effective potential energy between opposing unit statistical charges separated by 
a distance R is given by V{R) = — A;BTln(expz[0(O) — (f){R)]), where (f){R) denotes the 
phase variable of the dual XY-model. ■'^^'■'^^ Yet since the correlation function for this model 
is given by (expz[(^(0) — (piR)]) = (ro/-R)^/^exp(— i?/^^) near Tc, we have that the string- 
tension [defined by limi^^oo^(-R) = o'R] is cr{T) — /cbTc^~^(T) in this region. [This result 
agrees with Eq. (6).] On the other hand, the zero-temperature string-tension in the strong- 
coupling limit is given by cr(0) = ^gQ?iu}oa~^ , since a confining flux-tube corresponds to 
a smooth step in the DG-model for surface- roughening (7). In fact, it is evident that the 
step free-energy per length of this model is equivalent to the string-tension in general. The 
temperature-dependence of the former quantity is well known,^^'^^ and agrees with the 
results just cited. 

Last, the specific-heat per site for the interface roughening model (7) obtained from 
Monte-Carlo simulations^^ is reproduced in Fig. 1. It is qualitatively similar to the 
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weak-coupling result discussed previously, with the exception that (i) the strong-coupling 
specific- heat saturates to the classical value of above the transition because ksTc ^ 
hojQ, and that (ii) the energy-scale for activated behavior is on the order of a{0)a = 
0.69/cbTc. The latter observation confirms the fact that the elementary excitations in 
compact QED3 at low-temperature are strings of electric fiux. 

Anyon Superconductors and Spin-liquids. We now apply these results to the spin-1/2 
anyon superconductor,^^ which is known to be a large- A?" saddle-point of the 2D t — J model 
near half- filling. Consider then ideal spin-1/2 psuedo-fermions on the square-lattice with 
a statistical flux-tube of one flux-quantum attached to each. The corresponding long- 
wavelength limit CS lagrangian is given by 

L = ^|^(/za/za - c.c.)-(tV[^°+«]-«/,,/,+^,, + c.c.) + {27T)-h^,a^{i)a,{i)+ 

+ (Snyhoalii) - (87r//o)"'[Aia2(i) - A2ai(i)]'| (8) 

in the temporal gauge 03(2) = 0, where t denotes the hopping matrix-element for the 
psuedo-fermion (fia) and ai,{i) = hdai,{i)/dt, etc.. Also, eo and fio respectively denote 
the vacuum dielectric constant and the magnetic permiability for fluctuations a^(z) of the 
statistical gauge-field, -"^^'-^^ while Aq describes the statistical magnetic field due to the flux- 
tubes averaged over the entire lattice. After integrating out the psuedo-fermions, one 
finds that the cancellation of the resulting Hall conductance against that of the CS term 
leads to the action (1) for compact QED3 at zero-temperature, ^^'^^ with a zero-sound veloc- 
ity cq = /{eo + Sfy/'^ and a coupling-constant = {2n)~^{eo + ef)hu;o. 
Here, e/ ~ Aj^ and /ij^ ~ ^f^} psuedo-fermion contribution to the dielectric con- 
stant and the magnetic permiability, respectively, with Af and ^/ denoting the relevant 
(Hofstadter) energy-gap and magnetic length of the psuedo-fermions in the mean-field 
approximation. ^-"^ Now suppose that we extend this result above zero-temperature. Then 
no free statistical charge can exists below Tc because of confinement, implying that the 
Hall conductance generated by the psuedo-fermions is temperature-independent. Due to 
the exact cancellation of the latter against the CS-term, we (z) recover the assumption that 
the statistical gauge-field is described by vacuum QED3 below Tg, and find (ii) that the 
electromagnetic response is that of a superconductor in this regime. ^^'^^ Hence, Tc is in 
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fact the critical temperature for lattice anyon-superconductors. This is in sharp contrast 
to straight-forward extensions of RPA methods to non-zero temperature,^^ which find that 
the critical temperature vanishes due to a spurious temperature dependence acquired by 
the psuedo-fermion Hall conductance. Zero-temperature RPA calculations find,^^'^^ how- 
ever, that the London penetration length varies as \l oc (eo + efY^"^. By extension, A^^ 
should have an s-wave-like temperature dependence at low-temperature, while it should 
jump to zero at the transition. In conclusion, particle-hole excitations in lattice anyon 
superconductors are bound by a confining string.'^' The similarity between the duality 
of this transition to that of the 2D XY model^^'^^ and the duality of conventional su- 
perconductors in the presence of fluctuating magnetic fields to the three-dimensional XY 
model^^ is striking. Also, if we assume that ksT^ <^ A/, then the specific-heat should 
be dominated by the gauge- field contribution shown in Fig. 1. Such "rounding" of the 
specific- heat anomaly actually occurs in certain high-temperature superconductors.^^ 

In the case of the quantum-disordered (spin-liquid) phase of the 2D quantum antifer- 
romagnet, on the other hand, CPi calculations indicate that chiral spin-fiuctuations® are 
also described by the action (1) for vacuum QED3, where cq is identified here with the spin- 
wave velocity.^'^ The coupling-constant is given by Qq = 12ti IS. f /Tluq in such case, where 
A/ denotes the spin-gap.^ Hence, spinous are confined below Tc, where the elementary 
excitations are again strings of statistical electric fiux. Also, the existence of a universal- 
jump in the long- wavelength autocorrelation for the electric fields implies a corresponding 
jump A((5^)/jV = (27r)~^5rQ//?c in the fluctuation per (jV) site of the total number of 
skyrmions (quanta of chiral spin) Q-^ = {27r)~^ '^^€ijAiaj{x) at the transition.'^'^-'^ Last, 
the gauge-flelds will contribute a smooth bump anomaly to the specific heat below (see 
Fig. 1), just as in the previous case. 

To conclude, it has been demonstrated that the CD transition shown by compact 
QED3 is dual to the 2D CG transition in both the weak and the strong coupling limits. It 
has also been argued that the corresponding transition temperature coincides with the crit- 
ical temperature of ideal anyon superconductors and of spin-liquids on the square-lattice. 
The former constitutes a first step towards the thermodynamics of anyon superconductors, 
which are perhaps the paradigm for 2D superconductivity in general. 
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Figure Caption 



1. Shown are Monte-Carlo results obtained by Swendsen (ref. 25) for the specific heat of 
the DG model (7) on a 10 x 10 lattice. The specific-heat is plotted in units oiks, while 
the temperature is plotted in units of gohua/kB- The circles represent extrapolations 
based on the solid-on-solid model (ref. 25). 
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